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Ap in an SU(2)xU(l) chiral Yukawa model* 
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We report on a simulation in progress of an SU(2) xU(l) fermion-Higgs model with two non-degenerate fermion 
doublets. The Zaragoza proposal for chiral lattice fermions is used; naive fermions are studied for comparison. 
The aims of this simulation are a non-perturbative study of the decoupling of species doublers, of triviality bounds 
on the Higgs and fermion masses and of the phenomenological quantity A p. 



1. INTRODUCTION 

The study of fermion-Higgs models is moti- 
vated by the attempt to define chiral gauge theo- 
ries, and the Standard Model in particular, on the 
lattice. One assumes that corrections due to the 
electroweak gauge couplings can be incorporated 
perturbatively. The well-known obstacle to the 
inclusion of gauge fields is the species doubling of 
fermions, for which in the case of chiral theories 
no satisfactory solution is available [1]. 

In the Zaragoza proposal for chiral lattice 
fermions, the fermion interactions are constructed 
such that the species doublers are kept massless 
and non-interacting. The global chiral symmetry 
of the action is preserved (but, in the presence 
of gauge fields, gauge invariance is broken). This 
property is shared by the fermion action in the 
Roma I proposal, which can be shown to reduce 
to an action of the Zaragoza class [2]. 

Important issues to be investigated in chiral 
Yukawa models are the triviality of the Yukawa 
and Higgs couplings and the associated mass 
bounds. One would like to know in particular 
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if these remain within the perturbative domain, 
in a model with realistic fermion content. 

Here we report on a Monte Carlo study, in 
progress, of a fermion-Higgs model with two dou- 
blets of non- degenerate (Zaragoza) fermions and 
SU(2)xU(l) chiral symmetry. Apart from trivi- 
ality issues, we study the phenomenological quan- 
tity Ap, which depends on the fermion mass split- 
ting and whose leading contributions survive in 
the gaugeless limit of the Standard Model [3]. Its 
dependence on the fermion doublet multiplicity 
makes this quantity an interesting instrument to 
test the decoupling of species doublers. 

2. Ap 

The p parameter is defined as the relative 
strength of neutral and charged current interac- 
tions at zero momentum transfer, as measured 
in neutrino scattering. In the Standard model, 
these interactions are mediated by gauge boson 
exchange, and at tree level one has the relation 
p = M&/(cos 2 9 w Ml) = l. 

Quantum effects cause p to deviate from 1: 

Ap = 1-p- 1 . (1) 

The dominant contribution, from ("oblique") cor- 
rections to the W and Z propagators, is indepen- 
dent of the scattering process and is called the 
"universal part" of Ap. Additional contributions 
come from vertex corrections and box diagrams. 
They are small and process-dependent. Restrict- 
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ing ourselves to the universal part, p can be writ- 
ten in terms of the propagators of the W and Z 
bosons at zero momentum, 
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where Hw,z{p 2 ) are the transverse parts of the 
vacuum polarization tensors of the W and Z . 

The leading one-loop contribution is due to the 
top-bottom mass splitting. It takes the form [4] 
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if the bottom mass mj is neglected compared with 
the top quark mass m t = yv (with Gp 1 = \/2v 2 ). 
N c is the number of QCD colours. QCD correc- 
tions lower the right hand side by about 10 %. 

This result is a manifestation of "non- 
decoupling" [4]: the decoupling theorem [5], 
which says that observable effects of heavy parti- 
cles go as powers of their inverse mass, does not 
hold (at least, perturbatively) when the mass is 
generated by a Yukawa coupling. 

The experimental value of Ap constrains the 
top quark mass. The one-loop result (3) gives 
m t = 0(200) Gev, which is probably consistent 
with the perturbative treatment of the Yukawa 
coupling y. Alternatively, one could imagine Ap 
to saturate as a function of ra t beyond the per- 
turbative regime, allowing a larger ra t . For this 
scenario to be viable, the maximal fermion mass 
(Yukawa coupling) allowed by triviality would 
have to lie outside the perturbative regime. 

Apart from the leading result (3), there are 
"subleading" one-loop diagrams involving inter- 
nal gauge field lines, which are relatively sup- 
pressed by powers of g 2 /y 2 (which, by the way, is 
not a negligibly small number). In fact, at each 
loop level there is a leading term of zeroth order 
in g 2 , accompanied by subleading terms which 
vanish in the gaugeless limit g 2 — > 0. Hence all 
the leading contributions to Ap survive when the 
gauge fields are turned off. 

We shall now focus on this gaugeless limit. This 
has the great benefit that calculations can be 
performed within the fermion-Higgs sector of the 
Standard Model [3]. Gauge fields are only needed 
to act as sources for the currents whose two-point 



correlation function is calculated. From the Ward 
identities 
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where Zo t + are the wave function renormaliza- 
tion constants of the neutral and charged Nambu- 
Goldstone bosons (or "pions") tt°' + at zero mo- 
mentum [3]. We will use expression (5) to com- 
pute Ap non-perturbatively in a lattice Monte 
Carlo simulation. 

3. THE MODEL 

Consider the SU(2) x U(l) invariant fermion- 
Higgs model with action 

S = -|^Tr($+ + + h.c.) 
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Here $ is a (fixed-modulus) complex Higgs field, 
written as an SU(2) matrix, ipi^ are two (top- 
bottom) fermion doublets (two because we use 
the Hybrid Monte Carlo algorithm) and Y = 
U 







with y the top quark Yukawa cou- 



pling. The Yukawa coupling for the bottom 
quark(s) has been set to zero. The fermion field 
ipj^ix) in the interaction terms is defined as the 
average of the fields ipj at the 16 vertices of the 
hypercube at x, ensuring the decoupling of the 
fermion doublers from the path integral. As a 
consequence of the SU(2) x U(l) symmetry of the 
action (6), a calculation of the leading contribu- 
tion to Ap in this lattice model is guaranteed to 
produce the correct continuum result. This was 
demonstrated explicitly at one loop in Ref. [6]. 

For comparison we will also study the model 
with naive fermions, which amounts to replac- 
ing ipt 1 " 1 by tp in the action (6). In that case 
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Figure 1. Weak coupling region of the phase dia- 
gram for the model of Eq. (6). Both Monte Carlo 
data (diamonds) and mean-field results (dashed 
lines) are shown. 



the model describes 32 identical non-degenerate 
fermion doublets in the scaling region. 

4. TOWARDS RESULTS 

The relevant part of the phase diagram of the 
model with Zaragoza fermions is shown in Fig. 1. 
There are symmetric (PM) and broken (FM and 
AFM) phases separated by second order transi- 
tions. Monitoring of am t j av and am a / av (where 
v = 246 Gev sets the scale) in the scaling region of 
the FM phase will provide information on the up- 
per bounds on the top and Higgs masses imposed 
by triviality. By calculating Ap along the way, 
its dependence on ra t can be determined, from 
m t = up to the triviality "end point" . Here we 
shall limit ourselves to a discussion of the deter- 
mination of Ap and the problems involved. 

First of all, one would like to reproduce the 
weak-coupling result (3) numerically, in the scal- 
ing region. In Fig. 2 we compare the neu- 
tral pion propagators for naive and Zaragoza 
fermions. In both cases, the deviation from tree- 
level behaviour l/p 2 is very small (recall that 
Z cont = 2k Z\ a u). A salient difference is observed 
at p 2 = 4,8,12. As explained in Ref. [7], the 
"dips" in the naive case are enhancements of the 
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Figure 2. Monte Carlo results for the inverse 
lattice propagator of the neutral pion against 
p 2 = (1 — cosp^), in the small-t/ region. 



propagator due to the doublers. The fact that 
such dips are hardly observable in the model with 
Zaragoza fermions is a nice illustration of the ab- 
sence of doublers there, in the continuum limit. 

The charged pion propagators, not shown here, 
are virtually indistinguishable from the neutral 
ones, except at very small momenta. Hence a 
very careful analysis will be required to extract Zo 
and Z + from the propagator data with sufficient 
precision to determine Ap reliably. 

A poor man's method to obtain numbers for 
the Z's is to equate the data point at the smallest 
momentum value to p^^/Z. This is known [7] 
to overestimate the value of Z, because of the 
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curvature of the full inverse propagator at small 
momenta. Let us nevertheless try this for the 
data in Fig. 2, to get a first impression. It is 
interesting to compute the quantity 
Ap 
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using Eqs. (1,5), with dr = am t /av. N should 
come out as the multiplicity of fermion doublets, 
cf. Eq. (3) (32 for naive fermions, 2 for Zaragoza 
fermions). The poor man's method gives 224 ±26 
and 8 ±4, respectively. These numbers are far off, 
which is not surprising because Ap depends sensi- 
tively on the Z's. It is encouraging, though, that 
the ratio of these numbers is large, supporting 
the absence of doublers in the scaling region in 
the Zaragoza model. 

A much more accurate method, which has 
proven very successful in other fermion-Higgs 
models [7], is to fit the data with the propaga- 
tor in one-fermion-loop, finite- volume renormal- 
ized perturbation theory. As for now, let us just 
see what can be learnt from these one-loop self- 
energy corrections themselves. 

Fig. 3 is a plot of the (^-independent) quantity 

"Ap» n w+ ( P )-n w o(p) 
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y 2 /l6ir 2 p 2 (y 2 /I6ir 2 ) 

as a function of p 2 . n„.+,o(p) are the one-loop 
pion self-energy corrections due to one Zaragoza 
fermion doublet, calculated in a finite volume and 
with a finite lattice cutoff. This quantity is inter- 
esting because a Ap v (0) = Ap and "A/9"(p m ; n ) is 
essentially the poor man's estimate of Ap. 

This figure illustrates the difficulties in obtain- 
ing the correct value of Ap. First of all, the rapid 
decay (with a width ~ (am t ) 2 ) exposes the fail- 
ure of the poor man's method. Second, 'W'(O) 
is much larger than the value N = 1 expected in 
the continuum limit in infinite volume. This is a 
finite volume effect: doubling the lattice volume 
reduces "A/o"(p m ; n ) by a factor 2 (even though 
Pmi n is smaller in the larger volume!). Third, the 
finite cutoff leaves its marks: in infinite volume, 
"A"(0) = 0.82 for am t = 0.3, "A"(0) = 0.51 
for am t = 0.8. All these trends are similar when 
naive fermions are used. 

Problem #1 can presumably be avoided by us- 
ing the more sophisticated propagator fits men- 
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Figure 3. "N"(p) against p 2 , cf. Eq. (8). 



tioned earlier. Problem #2 may be overcome by 
extrapolations to infinite volume. Problem #3 is 
probably the most difficult to handle, since rela- 
tively large values of am t are inherent to the in- 
vestigation of triviality issues. On the other hand, 
a strong cutoff sensitivity (scaling violation) of 
Ap might actually teach us something about the 
cutoff in the Standard Model or even allow us to 
catch a glimpse of "the other side" . 
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